Nontopological Finite Temperature Induced Fermion Number 

I. J. R. AitchisonW and G. V. Dunne^ 1 ' 2 ) 
W Department of Theoretical Physics, Oxford University, 1 Keble Rd., Oxford 0X1 3NP, U.K. 
Research Centre for the Subatomic Structure of Matter, University of Adelaide, S.A. 5005, Australia 

We show that while the zero temperature induced fermion number in a chiral sigma model back- 
ground depends only on the asymptotic values of the chiral field, at finite temperature the induced 
fermion number depends also on the detailed shape of the chiral background. We resum the leading 
low temperature terms to all orders in the derivative expansion, producing a simple result that can 
be interpreted physically as the different effect of the chiral background on virtual pairs of the Dirac 
sea and on the real particles of the thermal plasma. By contrast, for a kink background, not of 
sigma model form, the finite T induced fermion number is temperature dependent but topological. 



The phenomenon of induced fermion number due to the interaction of fermions with topological backgrounds (e.g., 
solitons, vortices, monopoles, skyrmions) has many applications ranging from polymer physics to particle physics 
jj], ||, |J, ||, [| |7[ |i| . The original fractional fermion number result of Jackiw-Rebbi [Q has a deep connection with the 
existence of spinless charged excitations in polymers Q . The adiabatic analysis of Goldstone-Wilczek (j| in systems 
. without conjugation symmetry has important implications for bag models monopoles, and sigma models, which 
r ^ ' provide effective field theory descriptions of systems ranging from condensed matter, to AMO, to particle and nuclear 
physics Q. The induced fermion number is related to the spectral asymmetry of the relevant Dirac operator, and 
mathematical results concerning index theorems |^] relate the fermion number to asymptotic topological properties of 
the background. At finite temperature, the situation is less clear. In several examples |l0[ [ll], [l2|, |l3|, [uj, the fermion 
1 ■ number is known to be temperature dependent, but is still topological in the sense that the only dependence on the 
background field is through its asymptotic properties. In this Letter, we present a simple physical case for which 
i this is not true : in a 1 + 1 dim chiral sigma model, the finite temperature induced fermion number depends on the 
£SJ ' detailed structure of the background. This contradicts a previous analysis jl5) and claim |l6| that the finite T fermion 
£NJ , number is in general a topological quantity. We give a simple physical explanation of the origin of the nontopological 
dependence. Our analysis has been motivated in part by the results of |17| concerning the T dependence of anomalous 
, amplitudes in nuclear decays. 

Consider an abelian model in 1 + 1 dimensions with fermions interacting via scalar and pseudoscalar couplings to 
two bosonic fields </>i and 4>2- For the purposes of this paper </>i and </>2 can be considered as classical external fields. 
The Lagrangian is 



o 
o 



Ph| £ = i#V-V>(0i+*7502)^ (1) 

<D . 

' There are two especially interesting physical cases: 
, (i) kink case M : 

•rH , 

! <pi = m and 2 (±oo) = ±<f> 2 (2) 

(ii) sigma model case || : 

4>\+4>l= m 2 (3) 
In the sigma model case (||), the interaction term in the Lagrangian ([l]) can be expressed as 

m-i/i e 1150 tp = rwi/j (cos 9 + 175 sin 9) tp (4) 
At T = 0, both these cases have an induced topological current J' 1 =< ^^vb > given by || 

= _L e ^ dv6 + ... (5) 

Ztt 

where the angular field 9 is defined by 9 = arctan(02/</'i)- The dots in (^) refer to higher derivative terms, which 
are all of the form of a total derivative of 9 and its derivatives Q. Thus, in particular, the induced fermion number, 
N = JdxJ°, is 



(G) 
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where ±8 are the asymptotic values of 8(x) at x = ±00. The fermion number TV is topological as it depends only on 
8, not on the detailed shape of 8{x). The conjugation symmetric case of Jackiw and Rebbi jl[] is obtained by taking 



in -> in the kink case (0), in which case N — > ±| 



At nonzero temperature, the induced fermion number for a static background is Q 



2 J c 2ni \H - z J \ 2 

where j3 = 1/Tis the inverse temperature, and tr( „_ ) is the resolvent of the Dirac Hamiltonian H . The contour C is 
(— oo + ie, +oo + ie) and (+00 — ie, — 00 — ie). By considering static backgrounds we avoid the well-known complications 
of finite temperature calculations in non-static backgrounds |18| . The technical part of the calculation of the induced 
fermion number (Q) is the computation of the resolvent of H . Once this is done, the induced fermion number may be 
expressed as an integral representation, or as a sum by deforming the contour in (Q) around the simple poles of the 
tanh function. For static backgrounds 4>i ( x ) and 02 (x) in ([!]) , the Dirac Hamiltonian is 

H = -i 7 yV + 7Vi(aO+n°75<Mz) (8) 

where V = 4-, and we will work with the Dirac matrices 7 = 03, 7 1 = ia%, and 7 s = —o\. Also, note that only the 
even part (in terms of the argument z) of the resolvent tr( H 1 _ z ) contributes to the induced fermion number N in (Q). 
(This is most easily seen by deforming the contour around the poles of the tanh function.) 

Consider first the kink case in (Q). Then the even part of the resolvent can be computed exactly using a trace 
identity which is a special case of the Callias index theorem || [l(], |l^] (alternatively, it can be derived in a more 
elementary manner as an exact resummation of a SUSY derivative expansion p0|) : 



tr 



1 



H-z 



TO 

= tr I — — — — — I - tr 



-(V + 2 )(V - 0a) + to 2 - z 2 J V-(V - 2 )(V + 0a) + to 2 - z 2 

TO02 



(m 2 — z 2 ) \jm 2 + 02 — z 2 
Then the induced fermion number (0) for the kink case (||) is 



(9) 



AT = 2 - ( ^ ) sW £ 1 (10) 

* \ * / „ =o ((2 n+ l)2 + (!^)2)y / ( 2n+ l)2 cos 2 0+ (ir^)2 




FIG. 1: Plots of nN, where TV is the finite temperature fermion number ( |lo|) for the kink case (0), as a function of 9. These 
plots are for mf3/ir taking values 0.5, 1, and 10, as labelled. As T — -> 0, note that irN — > 6, as in (g). 
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where 9 = arctan(</>2 /m) . This result is consistent with previous analyses [jl(|, although these were much less explicit. 
The induced fermion number ( |To| ) is plotted in Fig. 1 as a function of 9 for various values of temperature. As T — > 0, 
this result reduces smoothly to the zero temperature result (^|) . Despite its complicated form, the nonzero temperature 
result ([l0|) is still topological as it only refers to the background through 6. 

In the sigma model case (|^) , the trace identity formulae (^) do not apply. Another approach is needed to evaluate 
the resolvent. One such approach is the derivative expansion pi] , in which we assume that the spatial derivatives of 
the background fields are small compared to the fermion mass scale m. In other words, the backgrounds (j>i(x) and 
4>2( x ) are assumed to be slowly varying on the scale of the fermion Compton wavelength. Returning to the general 
Hamiltonian (H), the derivative expansion can be obtained by separating H 2 as 



H 







cf)' 2 icj)' x 



(11) 



and then expanding tr y fJ 
yields: 



tr I (H + z) H 2_ 2 1 in powers of derivatives. A simple calculation to first order 



tr 



H-z 



dx ■ 



b 2 - z 2 ) 3 / 2 



(12) 



where the dots refer to terms with three or more derivatives. 

In the kink case (|2|), where <j>\ = m is constant , this first order calculation actually reproduces the exact trace 
identity result (^). But in the sigma model case (|J), where <j)\ + <j)\ = m 2 is a constant, the first order derivative 
expansion result (Il2|) implies that: 



tr 



1 



H-z 



2(to 2 - z 2 ) 3 ' 2 



dx9' 



(13) 



So, to first order in the derivative expansion, the induced fermion number for the sigma model case is 




' [(2n + l)2 + (^)2]3/2^ 



dx t 



(14) 



which is simply the zero temperature answer (|^) multiplied by a smooth function of T. As T — > 0, this prefactor 
reduces to so the full zero temperature result (||) is regained. But at finite temperature, the first order (in the 
derivative expansion) formula ( [14] ) for the sigma model case differs from the kink case formula ([lO]) , even though each 
of and @ reduces to (§) at T = 0. 

This raises the question of the higher order corrections to the derivative expansion (|T^). In the kink case (j^), 
there are no higher order corrections to the even part of the resolvent in (pd|). This is due to the special form of 
the Hamiltonian in the kink background, which leads to the first order formula (12) agreeing with the exact trace 
identity result (^|). There can, of course, be higher order corrections to the induced fermion number density, but these 
are all total (spatial) derivatives, and do not contribute to the integrated induced fermion number, even at nonzero 
temperature. 

But in the sigma model case (||), where the trace identity does not apply, the situation is very different. Going to 
the next order in the derivative expansion, we find 



tr 



H 



2(m 2 - z 2 f/ 2 



dx9' 



8(m 2 - z 2 f/ 2 
m 2 (Az 2 + m 2 ) 
16(m 2 - z 2 )V 2 



dx9"' 
dx(6 



l\3 



(15) 



where the dots refer to terms involving five or more derivatives. For a chiral background with 9{x) approaching its 
asymptotic values exponentially fast, the term J dxQ'" vanishes. But J dx (9 1 ) 3 does not vanish. Thus, the first order 
induced fermion number (Jhfl) acquires a third order correction: 



2 




[-4(2n+l) 2 + (^ 
' [(2n + l) 2 + (^)V /2 




/\3 



(16) 
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This is not just a function of the asymptotic value 9 of the chiral field 8(x); it also depends on the actual shape 
of 9{x). Thus, the induced fermion number is no longer topological. This contradicts [fl5| , where it is stated that 
the first order derivative expansion contribution (|l4|) is the full answer. However, the energy trace prefactor in ( |l6| ) 
vanishes at T — 0, so the nontopological third order contribution ([l6]) vanishes at T = 0. Thus, the nontopological 
nature of the finite temperature induced fermion number is still consistent (at this order) with the topological nature 
of the zero temperature induced fermion number (^) . 

We now turn to a physical explanation of why, in the sigma model case, the finite temperature induced charge 
is more sensitive to the background field than at zero temperature. Note first of all that the chiral background 
acts like a static but spatially inhomogeneous electric field, as can be seen by making a local chiral rotation [Q: 
ip — ► ip = e 1015 / 2 ^. In terms of these chirally rotated fields the Lagrangian (|l|), with interaction (|J), becomes 

_ _ _ Ql 

L = iippip — mipip — ip 7° — V" (17) 

(The chiral rotation leads to an anomalous Jacobian in the path integral, but this does not affect the induced fermion 
number.) Thus, the chiral field acts as an inhomogeneous Aq(x) — ^9'(x), leading to an electric field 

E{x) = \e"{x) (18) 

Given that 9{x) itself has a kink-like spatial profile, the electric field is such that it changes sign as a function of x, 
as shown in Fig. 2 (we choose 9' > 0). This electric field acts on the Dirac sea to polarize the vacuum by aligning 
the virtual vacuum dipoles of the Dirac sea, producing a localized build-up of charge near the kink center. But at 
nonzero temperature, the electric field also has an effect on the thermal plasma, as we show below. 

First, consider the full derivative expansion ( |l2] ) of the even part of the resolvent, at low but nonzero temperature. 
At fifth order, there are three independent terms involving 8""', 9"'(9') 2 , and (9 1 ) 5 . The 9'"" term vanishes when 
integrated over x, but the other two terms are generally nonzero. However, as T —* the (9 1 ) 5 term dominates the 
9"'(9') 2 term. Indeed, for low temperature, the dominant term with (21 — 1) derivatives in the derivative expansion 



( |l2| ) involves (9') 21 1 . Then, using the chirally rotated form (17) of the Lagrangian, the dominant term at (21 — l) th 
order is simply: 




with Euclidean p = (oj n , k) and u n = (2n + 1)ttT the Matsubara modes. 

At zero temperature, all these terms iV^ 2 ' -1 ) vanish, except for I = 1. This fact is not obvious; it involves highly 
nontrivial cancellations between terms in the expansion of the trace. But at nonzero temperature, all the terms in 




FIG. 2: For a kink-like chiral field 9(x) (dashed line), the corresponding electric field in (|1S|) has the form shown in the solid 
line, producing a vacuum polarization charge distribution localized near the kink center, roughly following the dotted line 8' . 
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( |l9| ) are non- vanishing. Moreover, they have a remarkably simple low temperature (T <C m) limit: 



2/-1 



N^ = ^j^-f^-^_L_ ijix ^ + ... (20) 

Thus, in the low temperature limits we can resum the entire derivative expansion, to obtain the induced fermion 
number in the sigma model case (F 



N=ldx^--J^ /^e-^sinhf^U... (21) 



2tt V 7T ./ \2T 

where the dots refer to subleading terms for T <^ m. 

Several features of this result ( pi] ) deserve comment. First, note that at zero temperature, only the first term 
survives, producing the familiar result (^) that the induced fermion number depends on the chiral field 9{x) only 
through its asymptotic value 9 = 9(oo). At zero temperature, one can invoke Lorentz invariance to constrain the form 
of higher order corrections to (^|) to be total derivatives but these arguments do not apply at finite temperature. 



We see this in (21): the temperature dependent corrections are not total derivatives of terms made from 9 and 
its derivatives. At nonzero temperature this shows clearly that the induced fermion number is nontopological - it 
depends also on the detailed shape of 9(x). Second, the resummed exponential factors e~ ( ™ =F6 ' / 2 )/ T in ( |2l| ) are 
consistent with the derivative expansion assumption that 9' <C m. Finally, the form of these exponential factors 
suggests an interpretation of the derivative expansion as an adiabatic change of the local Fermi level with a local 
chemical potential /i = — #'/2, which once again is only sensible in the derivative expansion regime where 9' <C m. 

To make this physical picture more precise, we can interpret the result ( pl| ) as follows. The first, topological, term 
refers to the induced charge coming from the polarization of the Dirac sea. This is temperature independent as the 
short-lived virtual "electron-positron dipoles" of the Dirac sea do not come to thermal equilibrium. The next term 



in (21) corresponds to the induced charge arising from the response of the real charges in the thermal plasma to the 
spatially inhomogeneous electric field (|l8|). Indeed, the linear response |2^] of the plasma at low temperature to such 
an electric field yields an induced fermion number density 

/dk 

where the static distribution function f(x,k) satisfies the Boltzmann equation 

v-^f(x,k) = -E(x)-^f(x,k) (23) 



where v — kj\Jk 2 + m 2 . Regarding fi — — ^9'(x) as a local chemical potential, ( p3[ ) is satisfied by local Fermi particle 
and antiparticle distribution functions 

f±(x, k) = ■ - (24) 

Inserting f = f + — /_ into (p2|), we obtain precisely the second, nontopological, term in ( pll) in the low temperature 
limit. 

At T=0, the fermion number may be defined as a sharp observable p^| ; but at T > 0, thermal fluctuations introduce 
an rms deviation. Thus, the finite T fermion number in (0j2l|) is a thermal expectation value (N), as in the monopole 
cases 0, |l3|, [l4|. We have estimated (N 2 ) — (A) 2 , in the derivative expansion regime, in an analogous manner to 
the computation presented here for (N). We find that the rms deviation vanishes at T=0, but at nonzero T can be 
significant compared to the thermal shift in (|2~l|). Details of this will be reported elsewhere. 

To conclude, we comment briefly on possible implications of these results for models in other dimensions for which 
there is an induced fermion number due to some nontrivial background. In 2 + 1 dimensions, fcrmions in a static 
magnetic background acquire an induced charge that is topological, expressed in terms of the net magnetic flux of the 
background. At finite temperature, the induced charge remains topological, but is multiplied by a smooth function 
of the temperature |pTf . In 3 + 1 dimensions, fermions in a static Dirac monopole background acquire an induced 
charge that is temperature dependent at finite T, but still only depends on the background through the total magnetic 
charge and the self-adjoint extension parameter Jl2|, 113]. A more interesting case is a static 't Hooft-Polyakov monopole 
background, which has a characteristic size scale. Consider, for example, the coupling 

Ant = i> (4 + 4> + hsm) i> (25) 
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where ip is an isodoublet fermion, is a static SU(2) 't Hooft-Polyakov monopole, and <f> is the corresponding static 
Higgs field. We have computed the finite temperature induced fermion number, using the 3 + 1 trace identity used 
in the zero temperature case [[24), and we find precisely the same expression (|l(|) as in the 1 + 1 kink case, with the 
identification 6 — arctan(0/m), where (f> is the asymptotic value of the magnitude \4>\ = y/cj) a (f) a of the Higgs field. 
Given that ([lC]) reduces to (||) at T = 0, this monopole result is consistent with the familiar zero temperature result 
H Q that the induced fermion number is proportional to 9 |^5| . The 3+1 dimensional analogue of the 1 + 1 sigma 
model case ||]4|) is the sigma model with coupling 

Ant = mip (tt + 1757? • f) 1p 

= m$ Q(<? + g ] ) + ^(g - 5^75) i> (26) 

where f are su(2) generators, the fields ttq and 7? are constrained by ttq + 7F 2 = 1, and the fields g in the second line 
are defined by g = ttq + in ■ f . At zero temperature, there is an induced topological charge density 

J ° = 24^" tr (d^digg-'djgg^dkg) (27) 

The corresponding zero temperature integrated charge is given by the winding number of the background field g at 
zero temperature. We conjecture that at finite temperature this induced charge will acquire additional nontopological 
contributions similar to those found here for the 1 + 1 sigma model case. 
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